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Renormalization of Massive Lattice Fermions 
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The renormalization of a general action for massive lattice fermions is discussed. The analysis applies for all 
m q a. Preliminary results for the self energy at one loop in perturbation theory are presented. 



1. INTRODUCTION 

This paper is a progress report of efforts to cal- 
culate the renormalization of fermion masses and 
bilinear currents in onedoop perturbation the- 
ory [1]. When these calculations are finished, 
they will permit a determination of heavy-quark 
masses, they will give one-loop predictions for the 
tuning of improvement parameters in the action, 
and they will give a one-loop guide for extrapo- 
lating matrix elements to the continuum. 

After specifying a general action in sect. 2, 
sect. 3 sketches an all-orders derivation for mass 
and wavefunction renormalization in terms of the 
fermion self energy. We have the self energy 
to one-loop for the simplest action, and present 
results with and without tadpole improvement. 
Sect. 4 discusses current renormalization. 

2. THE ACTION 

Consider the action S = So + Sb + Se + • • •, 

So = rnoJ^)-WcJ^)A^ 

+<;J${x)tD1>{x) (1) 
+ 1 / + lv)Do " (1 " l )D+]i>{x), 



Sb = —^cicbC J ip(x)S -B(x)ip(x) 



Se 



| ac sC / ip(x)a-E(x)ip(x). 



(2) 



(3) 



Special cases are the Wilson action [2], which sets 
r = (*=l,CB=CE = 0, and the Sheikholeslami- 



Wohlert action [3], which sets r = C, = 1, cb = ce- 
To remove lattice artifacts in general the parame- 
ters mo a, r, (*, cb and eg must all be adjusted [4]. 
In a non-relativistic setting, however, it is enough 
to adjust mod, eg, and eg [4]. 

3. THE SELF ENERGY 

The self energy is related to the 

momentum-space propagator by 



s- i (p) = s - i (p)-np), 



(4) 



where So(p) is the free propagator. In perturba- 
tion theory is the sum of all one-particle 
irreducible graphs. The po-Fourier transform 
C(t,p) = (27T)- 1 J dpo ey ot S{p) obeys 



C(t,p)=Z 2 (p)e- E P^Q- 



(5) 



where /. ( , is the energy of a one fermion state with 
momentum p, and Q is a Dirac matrix satisfy- 
ing (Q70) 2 = Qlo- The • • • denote multi-particle 
states, which are irrelevant here. 

The self energy has the decomposition 



(6) 



in Dirac matrices. For a Euclidean invariant cut- 
off C and Afi = A V// are functions of p 2 only. 
With the lattice cutoff, however, they are con- 
strained only by (hyper) cubic symmetry. For em- 
phasis it is convenient to write, say, C(po,p). 

To obtain an expressions for /. ( , and Z2 one 
carries out the po integration with the residue 
theorem. For arbitrary p the energy /. ( , is the 



2 



solution of the implicit equation 
1 + m a + \r(,p 2 a 2 — C = 



cosh Ea + (l- A )Vl - V 2 sinh Ea, 
for E. The abbreviation 



V 



2 E,-(C-^ 



) 2 sin 2 



fl - An) 2 sinh 2 Ea 



(7) 



(8) 



Here the self-energy functions A ll {po,p) and 
C'(po,p) are evaluated at po = i-E". The solution 
of eq. (7), E = E p , defines the (lattice-distorted) 
mass shell of the fermion. The residue is 



Z^(p) = (1 - A ) cosh E p a - CVl - V 2 



-(Vl-V 2 + A ) sinh E p c 

T,j Aj) sin 2 Pj a 

(1 — Aq) sinh E p a 



(9) 



The notation / = (ia)~ 1 (df/dpo). In eq. (9) the 
self-energy functions A ll {po,p) and C'(po,p) are 
evaluated on shell, i.e. po = /' / . ( , . 

The p dependence of Z'lijp) is an artifact of 
the lattice cutoff. An acceptable definition of the 
wavefunction renormalization constant is 

z - 1 = e M 1 a _ Aq cogh Mia + A Q sinn Mi a - C* ( 1 0) 

at p = 0, where Mi = Eo is the (all-orders) rest 
mass of the fermion. 

Eq. (7) at p = determines the rest mass via 



e Mia = l + m a + A sinh M x a - C (11) 

2 - 

(C-A,) 2 



and the dynamic mass M2 = {d 2 E p /dp 2 ) p l via 



M ia 



M 2 a 



1 - A ) sinh Mia 



d 2 r 

-— \ A (iE P) p) sinh E p a - C(iE p ,p) 



(12) 



A total pi-derivative includes an explicit part and 
an implicit part through the /. ( , dependence. In 
eqs. (11) and (12) the self-energy functions and 
derivatives are evaluated at po = iM\ and p = 0. 

For a massless fermion, M\ = M2 = 0. The 
bare mass that induces Mi = obeys 

m 0c a = C(0,O; m 0c a). 
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Figure 1. Plot comparing M[ a vs. Mo a (solid 
curve ) and M[ 1] a vs. Mo a (dashed curve). The 
static point is indicated by the box [5]. 



The third argument of C denotes the parametric 
dependence. It is useful to take care of this term 
once and for all, and write 



e M ia = 1 + Mq(1 + Ao s i n h Mia -C* 



(14) 



where Moa = m a — m 0c a = (2k) 1 — (2k c ) 1 , 
and C{iM\a, 0; m a) = C{iM\a, 0; m a) — m 0c a. 

We turn now to one-loop results for r = (* = 1, 
cb = ce = 0, with and without tadpole improve- 
ment. In perturbation theory the rest mass has 
an expansion 



Mia = log(l +M a) i^^'Mfo. 



(15) 



1=1 



In the tadpole improved version Moa = Moa/uo, 
where «o is a suitable (gauge invariant) average 
link. In applications both «o and k c would be 
taken from Monte Carlo calculations. Below we 
choose «o = (8k c ) _1 - Figure 1 shows the one-loop 
correction to the rest mass Mi. As expected, m{ 1] 
is significantly smaller than m[^ . 

For the dynamic mass it is better to define a 
renormalization factor via M2 = ttI'iZm, where 

Moa(l + Moa)(2 + M a) 
7712(1 ~ 2C 2 (1 + M a) + r(M a{2 + M a) ( ' 

is the tree-level expression for the dynamic mass, 
except that the linear mass divergence is ab- 
sorbed, order by order, into Moa. The tad- 
(13) pole improvement is M2 = tH'^Zm, where ui2 is 



3 



0.25 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l 




0.05 1 





_q hi i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i r 

5 10 15 20 25 30 35 40 45 50 

M2C1 

Figure 2. Plot comparing Z^j (solid curve) and 
Z^ vs. M2CI (dashed curve). 

given by the right-hand sid^of eq. (16), but with 
Mo a 1 — y Mo a. The factors Zm have series 

< Z M = l + jTg 2 l< zV\ (17) 
1=1 

Figure 2 shows the one-loop renormalization of 
the dynamic mass. Again, \Z^\ is significantly 

smaller than l-Z^I- 

To define the perturbative coefficients for the 
wavefunction renormalization constant, factor 
out e Mia . The tadpole improved constant is 
Z'i = uqZ'z- The perturbative series are 

^"zUl + E^? 1 . (18) 
1=1 

The one-loop coefficient has an infrared diver- 
gence, which can be regulated with a gluon mass 
A. Figure 3 plots the IR-finite = Z^ 1] + 

log(A 2 a 2 )/(67r 2 ). Once again, z\^ is significantly 
smaller than z^ . 

4. VERTEX CORRECTIONS 

A full vertex function takes the form 

V(p\q) = S(p)T(p\q)S(q). (19) 

In perturbation theory T is given by the sum of 
all truncated three-point diagrams. To put the 
external lines on shell, one Fourier transforms 
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Figure 3. Plot comparing z\ vs. Mo a (solid 
curve) and z^ vs. Mo a (dashed curve). The 
static point is indicated by the box [5]. 

in po and go- Poles arise precisely as in the self- 
energy derivation, so the on-shell truncated ver- 
tex function is T(iE p ,p\iE q , q). Consequently, 
when normalization conditions introduce a factor 
1 + moa at tree level, the all-orders generalization 
is e Mia , just as with Z^ 1 ■ 
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